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Uhlmann's mixed state geomètric phase [Rep. Math. Phys. 24, 229 (1986)] is analyzed in 
the case of a qubit affected by isotropic decoherence treated in the Markovian approximation. It 
is demonstrated that this phase decreases rapidly with increasing decoherence rate and that it 
is most fragile to weak decoherence for pure or nearly pure initial states. In the unitary case, 
we compare Uhlmann's geomètric phase for mixed states with that occurring in Standard Mach- 
Zehnder interferometry [Phys. Rev. Lett. 85, 2845 (2000)] and show that the latter is more robust 
to reduction in the length of the Bloch vector. We also describe how Uhlmann's geomètric phase in 
the present case could in principle be realized experimentally. 
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I. INTRODUCTION 

The geomètric phase, first discovered by Berry |lj for 
cyclic adiabatic quantal states and generalized to non- 
Abelian 0, nonadiabatic 0, and noncyclic evolu- 
tions, has recently been suggested and demonstrated 
experimentally :6| as a tool to achieve quantum compu- 
tation that is resilient to certain types of errors. Thcsc 
analyses have been mainly concerned with the adiabatic 
geomètric phase as this phase has the advantage that it 
depends only upon controllable external parameters \J\. 
However, from another point of view, adiabaticity is a se- 
rious drawback as it means that such gates operate slowly 
compared to the corresponding dynamical time scale and 
this makes them potentially vulnerable to loss of coher- 
ence. 

This issue has recently been addressed by Fonseca 
Romero et al. |8j, who have analyzed the behavior of 
the geomètric phase in the presence of decoherence in 
the qubit (two-level) case and have shown that the con- 
tribution from the adiabatic geomètric phase to the pre- 
cession of the Bloch vector in a slowly rotating magnètic 
field is robust to certain kinds of anistropic weak damp- 
ing. Nazir et al. have analyzed implementations of 
geomètric quantum computation in the presence of de- 
coherence and have demonstrated a path-dependent sen- 
sitivity to anisotropic noise and loss of cntanglement for 
such gates. 

A different perspective on the effect of decoherence 
arises when considering geomètric phases for mixed 
states. Uhlmann (see also Refs. [lïllrllïllïl]) was 
probably first to address this issue by lifting the density 
operator to a pure state by attaching an ancilla and let- 
ting the combined system be parallel transported along a 
specific purified path. More rec ently , Sjòqvist et al. [Iq ] 
(see also Refs. EH E3, EH ES E3, IHf ) have discovered an- 
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other mixed state geomètric phase restricted to unitarily 
evolving density operators in the experimental context 
of one-particle interferom etry . Thegeometric phases for 
mixed states proposed in [lOj and [l5j are generically dif- 
ferent in the unitary case and match only under yery 
special conditions such as in the limit of pure states [23 • 
In this work, we analyze the behavior of Uhlmann's ge- 
omètric phase in the presence of decoherence and propose 
an experimental realization of this phase. Specifically, 
we consider this geomètric phase for a qubit affected by 
isotropic depolarizing decoherence treated in the Marko- 
vian approximation. The Markovian treatment of the 
depolarization channel makes it possible to continuously 
monitor the mixed state geomètric phase for the Bloch 
vector in its motion towards the origin inside the Bloch 
sphere. We hope that the present work would provide 
insights into Uhlmann's mixed state geomètric phase, ex- 
tcnding the purely mathematical treatments of this phase 
presented in the literature so far to an explicit physical 
situation. 

In the next section, we describe Uhlmann's approach 
to the mixed state geomètric phase. We apply the general 
formalism to the qubit case and derive the corresponding 
parallel transport equation first obtained by Hübner |13| . 
The unitary and Markovian approaches to the depolar- 
ization channel are described in Sec. III. Uhlmann's ge- 
omètric phase for a qubit evolution in the depolarization 
channel is computed in Sec. IV. In the idealized uni- 
tary pure state case, this evolution describes an isosceles 
spherical triangle on the Bloch sphere. We also describe 
how this geomètric phase could in principle be realized 
experimentally. The paper ends with the conclusions. 

II. UHLMANN'S GEOMÈTRIC PHASE 

The key idea in Uhlmann's ^3 approach to the mixed 
state geomètric phase is to lift the system's density oper- 
ator acting on the Hilbert space H to an extended Hilbcrt 
space 

n ext = n®n a (i) 
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by attaching an ancilla a. This extension has the prop- 
erty that every unit vector (purification) \tp) in the ex- 
tended space may be reduced to a density operator in the 
space of operators acting on TL as 



P = Tr|^|, 



(2) 



where Tr is partial trace over the ancilla. The lift of 
p is defined in terms of Hilbert-Schmidt operators W : 
TL a — > H such that p = WW^ . Such a lift is not uniquc 
as W i— > W — WU generates the same p for any choice of 
unitary U, i.e. we have a gauge freedom in the choice of 
U. Thus, we have an infinite number of ways to perform 
a purification of p. Our concern is now to pick out a class 
of exceptional purifications that defines a natural notion 
of geomètric phase. 

First, let us introduce the inner product between any 
pair Wi, W 2 of Hilbert-Schmidt operators as 



(W 1 ,W 2 ) =Tr[w}W 2 }. 
In terms of this we may define the quantity 



(3) 



v = ÇWi,W m ){W m ,W m -i) ■ ■ ■ {W 2 ,W X ) = £ (W u W m ) 

(4) 

corresponding to the ordered set II = pi,...,p rn of den- 
sity operators. Now, exceptional choices are the ones 
that maximize |£|, or, equivalently, that maximizc cach 
\(W j+u Wj)\. 

By making use of W — y/pV, V being unitary, and the 
polar decomposition y/p j+1 */p. = \ yfp~ 3 -\U j+ i ,j we 
obtain a necessary and sufficient condition for |£| maxi- 
mal when 



\(W J+1) W j )\ = Tt 



VP j+1 PjVp J+1 



(5) 



It follows that v — £ (Wi, W m ) is invariant under the re- 
maining gauge freedom Wj >— > ejUWj, \ej\ — 1 and U a 
fixed unitarity. Hence, with W being exceptional, v de- 
pends only upon the ordered set II and we define the cor- 
responding Uhlmann phase as 9 [II] = argz^. This phase 
is real-valued, gauge invariant, and, in the m — > 00 
limit, independent of the subdivision of the path. These 
properties make </> 9 [II] a natural definition of the mixed 
state geomètric phase associated with IT. Furthermore, 
by introducing <j)j+\j = arg(IUj + i, Wj) we may write 



í>3,2 + 1 



■tag(Wi,W„ 



(6) 



This may be simplified by choosing a particular class of 
gauge that corresponds to parallel lift defined by requir- 



mg 



W] +1 Wj > 0, Vi, 



(7) 



which implies (pj+ij — for each intermediate step be- 
tween 1 and m in Eq. ©. Note that the main condition, 



i.e. that |£| is maximal, is still true when Eq. holds. 
For this important class of gauge Uhlmann's mixed state 
geomètric phase reads 



FI] =arg(Wi,W„ 



(8) 



Let us now consider the parallel transport condition in 
the particular case of a qubit (two-level system). Any 
qubit state can be written asp=|(l + r-cr), where r 
is the Bloch vector, the length of which being less than 
unity for mixed states, and cr = (a x , <J y ,a z ) are the Stan- 
dard Pauli matrices. The parallelity lift condition Eq. 
(T7I) for any pair pj and Pj+i of such mixed qubit states 
|23| is equivalent to 



'j+r 



PjPj+WPj 



(9) 



By writing yfp] = ciq + a • cr and y/pj+ï — Òq + b ■ cr, and 
using features of 2 x 2 matrices we obtain 13] 

t = &oa + b.a + ? (bxa).cr 

3 V(Mo + b-a) 2 + |bxa| 2, 

For a continuous path II : t € [0, r] 1— > r(í) of qubit states 
we need to consider the infinitesimal version of Eq. 11(J|I . 
Let Vj = V, Vj+i = V + dV, y/pj+ï = ^fp] + d^/pj and 
use that a 2 , + |a| 2 = ^ we obtain to first order in da the 
diffcrential equation |13| 

dVVi = 2i(d& x a) • cr. (11) 

Formally we may write the solution of Eq. (|ll|l as 



V = V exp [li J(dax a) • cr 



Vo, 



(12) 



where V stands for path ordering and Vq is the initial 
unitarity. Evaluating this path ordered expression and 
inserting into Eq. JSJ) yields Uhlmann's mixed state geo- 
mètric phase for any qubit state. 

In the qubit case, it has been argued 0. Ull that the 
interior of the Bloch sphere is curved. This result is es- 
sentially captured by the mixed state line element |25j| 



ds 2 = 



dr 2 



1 



r 2 dn ■ dn 



(13) 



with r = rn, which shows that as one moves away from 
the center at r = 0, the circumference of the 2-sphere 
defined by each fixed r > grows more slowly than the 
distance from the center, due to the factor 1/(1 — r 2 ) 
in front of dr 2 . Uhlmann's geomètric phase for a qubit 
is the holonomy that naturally measures this curvature 
inside the Bloch sphere, just as the Standard pure state 
geomètric phase is the holonomy that measures the cur- 
vature of the Bloch sphere. In view of this and due to 
the non-Euclidean behavior of the radius and circumfer- 
ence when moving away from the origin inside the Bloch 
sphere, there is no reason to expect that the Uhlmann 
phase should have any direct relation to the sòlid angle 
enclosed by the Bloch vector in three dimensional Eu- 
clidean space. The result of the calculation in Sec. IV 
below for a qubit in the depolarization channel may be 
regarded as an illustration of this intuitive reasoning. 
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III. DEPOLARIZATION CHANNEL 

In the depolarization channel, the environment induces 
isotropic errors in the qubit state. This may be repre- 
sented by taking the three errors 

(i) |V») ^a x \tp), bit fl-ip, 

(ii) |V>) i-> C7- J »|V>), phase flip, 

(iii) \tp) i— > (j y \ip), both bit and phase flip, 

to be equally likely, each occurring with probability p/3. 
A unitary representation of the channel, using a minimal 
set of ancilla states {|0 a ), . . . , |3„)}, is given by 



U : p®\Ò a )(Ò a \ 



- \2 a ) ® <Ty ~ 



I 



|3.) 



xRprf y/l-pl® (Ó,| 

+ yf (ïa| + (2a| +cr z <8> (3«| 



(14) 



where R is some unitary operation acting only on the 
qubit. This may be lifted into the pure state evolution 
by adding another ancilla system b such that the initial 
state reads 



|6«><ò«| = 




and extending U to U €5 J&. 

Similarly, under certain restrictions one may model the 
map of the system's density operator as a nonunitary 
Markovian evolution described by the Lindblad equation 
(h = 1 from now on) 



P : 



(16) 

Each term L^pLj^ represents one of the possible er- 
rors (quantum jumps), while the sum over — ^L^L^p — 
TjpL^L^ is needed to satisfy the normalization condition. 
In this framework, the depolarization channel may be 
represented by the Lindblad operators 



where T is the time-independent decoherence rate. In- 
scrting these L M 's into Eq. I|l(jf> . we obtain the Lindblad 
equation for the depolarization channel as 



i[p,H] 



2L 



(18) 



Assuming the Hamiltonian H — ^lü(J z , corresponding to 
the unitarity R = exp ( — ^wícr z ), and the initial con- 
dition r(0) = ro(sin#, 0, cos#), the solution of Eq. Ijl8(l 
reads 



r(t) = r Q 



1 (sin 9 cos cot, sin 9 sinwí, cos 9). (19) 



Thus, the Bloch vector precesses uniformly around the 
z axis and its length decreases isotropically. The ef- 
fects of the unitary and Markovian treatments of the 
depolarizing channel can be formally related at time t 
as p(í) = f(l - r(i)/ro) = f (1 - e~ 4rt / 3 ). 



IV. UHLMANN'S GEOMÈTRIC PHASE IN THE 
DEPOLARIZATION CHANNEL 

A. Theoretical analysis 

To evaluate Eq. I|12l) analytically in the depolarization 
channel, we may choose a path where (da x a) • <r is time- 
dependent in such a way that it commutes at different 
times within a set of time intervals. Such a path is 

4r + 

A — > B : r(t) = r e~~ (sinwí, 0, coswí), 
< t < tt/(2w), 

4r j. 

B — » C : r(í) = r (sinwí, — coswí, 0), 

tt/(2lü) < t < (íp + tt/2)lü, 

4r . 

C^D : r(t) = ro e~~ (cos(/?sin[(jjí — ip], 
sin ip sin [u>t — ip], — cos [u>t — (p] ) , 
[tp + tt/2)/uj < t < (ep + n)/uj, (20) 

which is shown in Fig. 1 in the particular case where 
ip = n/2. 

The idealized pure state evolution is obtained when 
r = and ro = 1; it is represented by the sòlid line 
in Fig. 1 that defines the geomètric phase —w/A. In the 
general case, the third rotation is taken around the direc- 
tion m = (sin (p, — cos <p, 0) and the pure state geomètric 
phase becomes —<p/2. 

To compute Vd and thereby Uhlmann's mixed state 
geomètric phase for any fixed T, we first note that with 
a(t) — \&(t)\ it follows that da x a = a 2 (t)dn x n and for 
the above path we have 



H= 1,2,3, 



(17) 



A- 


-> B 


dn x 


n = 


— ujdt y, 


B - 


-> C 


dn x 


n = 


— íüdt z, 


C- 


4 D 


dn x 


n = 


—u>dt m 



(21) 
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FIG. 1: Path with istropically decreasing length of the Bloch 
vector. The sòlid line represents the unitary pure state case 
and defines a spherical triangle enclosing the sòlid angle 7r/2. 
The dashed line represents the case with nonvanishing deco- 
herence. 



where y, z, and m are unit vectors along the three rota- 
tion axes. Substituting this into Eq. H12[> yields 



r{ip+TT)/u) 

u = exp I — 2iuj / a 2 (t)dtm·cr 

^ J(<p+ir/2)/u> 

x exp ( — 2íüj / a2 (t) dt a z 

^ Jtt/{2lü) 

/.7T/(2w) 

x exp ( — 2iw I 

e -i X ™-<T e -ÍK<7 Se -inv y y Aí (22) 



a 2 (t)dt <j v )V a 



where we have used that Vq — Va- Now, from 
(y/p) 2 = p and y/p = a,Q + an ■ cr, we obtain a 2 {t) = 

^1 — y/l — r 2 (í) J /4 so that the angles fj,, k, and \ have 
the genèric form 

2uj / a 2 (í)dí = — 



1 — r\e 3 * 



(23) 



with ífc = 0, 7r/ (2a;), [y> + k/2]/uj and A</? = 7r/2, 7r/2, 
respectively. The integral can be evaluated by using 



+ ln 



1 — r^e 3 1 



= 2 



1 — a/1 — r 2 e 3 ' 



4a 2 (t) - 1 + ln 



a (í) 



a(t) 
(24) 



with the concomitant integration limits for /j, k, and x- 
Let us further introducc 



a (0) + a(0)a z = a + (3a z , 



PA 
~P~D 



«0 



(25) 



In terms of Eqs. (|22|) and (|25|l the geomètric phase for 
the path A — > £> — > C — » D becomes 

ç6 g = argTr[ v /pJ v /pB'í/ D V^] 

«77 + 



= — arctan 



av + /3tj 
sin k + sin[</5 — k] tan x tan /i 
cos k + cos[ip — k] tan \ tan /j 



(26) 



To further analyze this general analytic result, let us 
consider the following important special cases: 

1. Pure and unitary case: Here, r(t) = 1, V t > 0, 
which implies that a = (3 = v = T] = ^, p — 
X = tt/4, and k = (p/2 yielding cf> g = —(p/2. This is 
the expected pure state geomètric phase for a qubit 
enclosing the sòlid angle ip on the Bloch sphere. 

2. Mixed and unitary case: Here, T — so that r(t) = 
ro ^ 1, V í > 0, which implies that pjj = pa, 



X = tt(1 - V 1 -r'è)/ 4 = C, and k = ip(l 



y/1 — Tq)/2. Inserting this into Eq. l]2(:i[) we obtain 



arctan 



/ sin k + sin[(£> — k] tan 2 £ 
Vcos k + cos[ip — k] tan 2 £ 



(27) 



This case has also been analyzed in the context of 
Mach-Zehnder interferometry in [15(, where it was 
shown that the interference pattern for ro > 
is shifted by the mixed state geomètric phase 
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arctan 



ro tan 
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(28) 



where the sòlid angle cncloscd by the Bloch vector 
is ip. The mixed state geomètric phases in Eqs. Q27JI 
and H28JI match only if tan£ = 1 or K — (f/2 = nir, 
n integer. The former case corresponds to pure 
states and holds for any ep. As < k < (p/2 for 
ro > 0, the latter case has solution only for n = 
so that k — (p/2, which again corresponds to pure 
states. Thus, for the present evolution, (j> g and j g 
only match for pure states. In fact, one may show 
that 



sin k + sin[(p — k] tan 2 ( 



cos k + cos[(p — k] tan Ç 



< 



tan 



9 



,1 < l7 ff l 



(29) 
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with equality in the limit of puré states. This sug- 
gests that the mixed state geomètric phase pro- 
posed in 01 is more robust to reduction of the 
length of the Bloch vector than that proposed by 
Uhlmann. 

3. Pure initial state: This case is characterized by a = 
(3 = i, which yields 



arctan 



sin k + sva[íf — k] tan \ tan /i 



cos k + cos[ip — k\ tan \ tan /i 



(30) 



Thus, from Eq. ^ it follows that \<p g \ < (p/2 
for pure initial states. This reduction of the geo- 
mètric phase value as a function of the decoherence 
efficiency parameter Y/lo is illustrated in Fig. 2 
in the case where ip = n/2. We see that the geo- 
mètric phase decreases rapidly with increasing de- 
coherence rate. In particular, we note that when 
the decoherence rate is as large as the precession 
time-scale (i.e. Y/u = 1), the geomètric phase has 
practically vanished psj ]. 




2w 



\t)dt = 



-O 



(32) 



Comparison of Eqs. (|31|) and (|32|l shows that 
Uhlmann's mixed state geomètric phase in the case 
of weak isotropic decoherence is most fragile for 
pure or nearly pure initial states. 



B. Proposed experiment 

To test the above predictions experimentally requires 
control over the dynamics of the ancilla [23] ■ This is very 
difficult in the Markovian case, where the ancilla typi- 
cally has a complicated structure with many degrees of 
freedom so as to quickly forget the information acquired 
from the qubit. However, one could imagino a few-qubit 
implementation adapted to the unitary representation of 
the depolarization channel, in case of which Uhlmann's 
geomètric phase could be realized interferometrically. 

The general idea behind such a realization is based 
upon purification lift of the path 1 1— ► p(t) by adding an- 
cilla systems in such a way that partial trace over these 
ancillas is p(t). Such a measurement scheme may be con- 
structed as follows: 

(i) Represent the ancilla states /x = 0, . . . 3, with 
a pair of qubits so that |0 a ) = |0 a ) <£> |0 a ), |l a ) = 
|0„) ® |2„) = |1„> <8 |0 o ), and |3„) = |1„) ® |1„) 
and write the initial state Eq. I|15|l as 



\*a) 



'l + r A 



|0) ® |0 b ) 



FIG. 2: Uhlmann's mixed state geomètric phase (vertical 
axis) as a function of the decoherence efficiency parameter 
r/w (horizontal axis), V and w being the decoherence rate 
and precession frequency, respectively. The precession angle 
around the z axis is restricted to tt/2. 



4. Small decoherence rate Y/u> <C 1- First, in the case 
where ro = 1 we may expand the right-hand side 
of Eq. (|24|l to lowest order in the small quantity 

' S 3 t yielding 



v/T 



2w 



t k +Aip/uj 



O 



(31) 



On the other hand, when < rp 7^ 1 



and r^(8r/3)í max « 



1 ~ r 2 , 



one may take 



1/1 - rl + 0[Yt] so that 



'l-r A 



1)® |U) ® |0 a )<8) |0„), (33) 



where r A — r(0), and = |0) according to Fig. 
1. 

(ii) Take through the depolarization channel 

around the path of Fig. 1 yielding the final state 
in Schmidt form 



l + r D 



|0) ® \D) 



1 - l'D 



|1) ® \D A 



(34) 



with rjj — r([(p + vr]/w) and (D\D^) = 0, where 
\D), \ D ± ) ^TLb®TL a are normalized and read ex- 
plicitly 



6 



\D) = 




y/l-p D \Ob) <8 |0„) ® |0„) 



+ \/ í y|O b )® |1„><8>|1«) 




6 )(8|0 o )(8|l B ) 



+M/^|O b )® |la)® |O a ) 



Vl-Pl>|l6>® |0 o )O|0a) 



Ç 1 l„; l„; 



where p D = |(1 - r D /r A ). 



(35) 



(iii) Perform a unitary transformation under which [2í 



|£>) |0 6 ) ® |0„) <8 |0 o ) 
|^> -> \h) ® |0„) <8 |0 o ) 



(36) 



so as to obtain the state 



'l + r D 



|0) |0fc> 



1 - 7\D 



1) ® II 




0„)®|0 o ). (37) 



(iv) Expose \^d) to VoVl_®/6®/ a and let the resulting 
state interfere with e 1$^), 5 being a variable U(l) 
shift. The intensity X reads 



l^(l+r A )(l + r D )(0\V D VÍ\0) 

+ \y/0--r A )(l-r D )(l\V D Vl\l) 
Tr[y/pIy/pEV D vl]. (39) 



Thus, Uhlmann's geomètric phase for the path in 
Fig. 1 is realized as the Pancharatnam relative 
phase fio] that shifts the interference oscillations 
obtained by applying a variable U(l) shift to one 
of the interfering states. Such an experiment could 
in principle be implemented using, e.g., ion traps 
or NMR techniques, by letting three qubits act as 
the ancilla systems a and b for a fourth qubit in the 
reduced mixed state p. 



V. CONCLUSIONS 

We have computed Uhlmann's mixed state geomètric 
phase for a qubit affected by the depolarization channel 
and have described how this phase could in principle be 
realized experimentally. A ràpid decrease of this phase 
with increasing decohcrence rate has been demonstrated. 
For weak decoherence we have found that Uhlmann's geo- 
mètric phase is most fragile for pure or nearly pure states. 
In the unitary case, we have also demonstrated that the 
mixed state geomètric phase proposed by Uhlmann seems 
to be more sensitive to reduction of the l eng th of the 
Bloch vector than that proposed in Ref. [Ï5j. In this 
context, it would be interesting to compare the results of 
the present analysis with the extension of the mixed state 
geomètric phase in to the decoherence case [3l]. We 
hope that this work may lead to further studies of the 
Uhlmann phase for various quantum channels as well as 
to experiments of mixed state geomètric phases in the 
presence of decoherence. 



T oc 



D 



= 2 + 29? 
where 



(^A\V D V\ ® l6®/ a |*D) 



{^A\V D Vl®I b ®I a \^ D )e-' l& } 1 (38) 
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